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Abstract 

Simplified solutions to determine binary elements by astrometry were obtained 
in terms of elementary functions (Asada et al. 2004), and therefore require neither 
iterative nor numerical methods. In the framework of the simplified solution, this 
. . paper discusses the remaining two parameters of the time of periastron passage and 

CL|. the longitude of ascending node in order to complete the solution. We thus clarify a 

O ' difference between the simplified solution and other analytical methods. 

■ Key words: astrometry — celestial mechanics — stars: binaries: general 

^ ' 1. Introduction 
(N ■ 

Recently, we developed a formulation for determining binary elements with astrometric 
. observations. The simplified solution is written in terms of elementary functions, and therefore 
requires neither iterative nor numerical methods (Asada et al. 2004). This solution has been 
O ■ generalized to a binary system in open (hyperbolic or parabolic) orbits as well as closed (elliptic) 
ones (Asada 2007). An extension to observational data has been also discussed (Asada et al. 
^ ] 2007). The solution gives an explicit form of binary elements such as the eccentric anomaly 
5^ I and the major axis of elliptic orbits. Oyama et al. (2008) made an attempt to use this solution 
for discussing some uncertainty in binary elements because of large scatter of their data points, 
when they measured proper motions of maser sources in the galactic center with VERA. 

On the other hand, the remaining parameters of the time of periastron passage and the 
longitude of ascending node are not discussed in the simplified solution. Hence, the solution 
is rather simplified. However, these parameters are needed to make a comparison between the 
simplified solution and conventional ones. In addition, the lack of information on the remaining 
parameters apparently suggests a certain incompleteness of the simplified solution. In this brief 
article, therefore, we shall derive, in the framework of the simplified solution, both the time of 
periastron passage and the longitude of ascending node in order to complete the solution. 
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Astrometry plays a fundamental role in astronomy through providing useful star cata- 
logs based on precise measurements of the positions and movements of stars and other celestial 
bodies. For instance, astrometric observations provide an useful method of determining mass of 
various unseen celestial objects currently such as a massive black hole (Miyoshi et al. 1995), an 
extra-solar planet (Benedict et al. 2002) and two new satellites of Pluto (Weaver et al. 2006). 
Astrometry of Sharpless 269 with VERA detects a trigonometric parallax corresponding to 
a distance of 5.28 kpc, which is the smallest parallax ever measured, and puts the strongest 
constraint on the flatness of outer rotation curve (Honma et al. 2007). Accordingly, astrom- 
etry has attracted renewed interests, since the Hipparcos mission successfully provided us the 
precise catalog at the level of a milliarcsec. In fact, there exist several projects of space-borne 
astrometry aiming at a accuracy of a few microarcseconds, such as SIM^ (Shao 2004), GAIA^ 
(Mignard 2004, Ferryman 2004) and JASMINE^ (Gouda et al. 2007). 

In this paper, we focus on an astrometric binary, for which only one of the component 
stars can be visually observed but the other cannot, like a black hole or a very dim star. In this 
case, it is impossible to directly measure the relative vector connecting the two objects, because 
the secondary is not directly observed. The position of the star is repeatedly measured relative 
to reference stars or quasars. On the other hand, the orbit determination of resolved double 
stars (visual binaries), which are a system of two visible stars, was solved first by Savary in 
1827 and by many authors including Kowalsky, Thiele and Innes (Binnendijk 1960, Aitken 1964 
for a review on earlier works; for the state-of-the-art techniques, e.g., Eichhorn and Xu 1990, 
Catovic and Olevic 1992, Olevic and Cvetkovic 2004). The relative vector from the primary 
star to the secondary has an elliptic motion with a focus at the primary. This relative vector 
is observable only for resolved double stars. 

In conventional methods of orbit determination, the time of periastron passage is one of 
important parameters because it enters the Kepler's equation as 

t = to + ^{E-eKsmE), (1) 

where Iq, T, ck and E denote the time of periastron passage, orbital period, eccentricity and 
eccentric anomaly, respectively (e.g., Danby 1988, Roy 1988, Murray and Dermott 1999, Beutler 
2004). The simplified solution does not use the Kepler's equation in order to avoid treating 
such a transcendental equation. 

This paper is organized as follows. Our notation in the simplified solution will be sum- 
marized in § 2. The time of periastron passage in the simplified solution will be derived in § 3. 
The longitude of ascending node will be obtained in § 4. 

^ http://sim.jpl.nasa.gov/ 

^ hthttp://www.rssd.esa.iiit/index.php?project=GAIA&page=index 
^ http://www.jasmine-galaxy.org/ 
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2. Simplified solution: Our notation 



Our notation in the simplified solution is briefly summarized as follows. We neglect 
motions of the observer and the common center in our galaxy. Namely, we take account only of 
the Keplerian motion of a star around the common center of mass of a binary system. Let us 
define {x,y) as the Cartesian coordinates on a celestial sphere, in such a way that the apparent 
(observed) ellipse on the celestial sphere can be expressed in the standard form as 



a 



where a > b. The eccentricity e is yl — b'^/a'^. This eccentricity may be different from eK, 
the eccentricity of the actual elliptic orbit, because of the inclination of the orbital plane with 
respect to the line of our sight. The star is located at Pj = {xj,yj) on the celestial sphere at 
the time of for j = 1 , ■ • ■ , n. 

We use a fact that the law of constant-areal velocity still holds, even after a Keplerian 
orbit is projected onto the celestial sphere. Here, the area is swept by the line interval between 
the star and the projected common center of mass but not a focus of the apparent ellipse (See 
Fig. 1). This fact is expressed as 
S_ ^ S{k,j) 

T nk,jy 

where S{k,j) and S denote the area swept during the time interval, T{k,j) = —tj for tk > tj, 
and the total area of the apparent ellipse nab, respectively. The swept area is expressed as 
(Asada et al. 2004, Asada 2007) 
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. . . , y. 



S{k,j) = -ab Uk — Uj -{smuk — sinu^) + ^(cosn^ — cosuj) . (4) 



a' " " b 



The eccentric anomaly in the apparent ellipse is given by Uj = aicta.n{ayj /bxj) . 

The orbital elements can be expressed explicitly as elementary functions of the locations 
of four observed points and their time intervals (Asada et al. 2004). Let us take four observed 
points Pi, P2, P3 and P4 for ti < ^2 < ^3 < ^4- The location {xe,ye) of the projected common 
center is given by 

F1G2 — G1F2 . . 

Xe = — Q ^ ^ ^ ^ , (5) 

-C/l-r2 — -ri£/2 

G1E2-E1G2 

= — w~F^' 

r^ir2 — riJ^2 

where Ej, Fj and Gj are elementary functions of T(j + 2, j + 1), T(j + 1, j) and Uk for k = j, 
j + 1, J + 2. The eccentric anomaly in the actual ellipse (on the orbital plane) is denoted as E 
(See Eq. (1)). 

Given a, b, Xe and ye-, we can analytically determine the parameters ck, i, clk and u as 
(Asada et al. 2004) 
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cosz = -(e-v/e-4), (8) 



where 



1 + cos^ i ' 

008207=^5 2—, (10) 

^=;;^V'^Hl/e^ (11) 



+ /i_e2 

3. Time of periastron passage 

In order to determine ax and for an actual ellipse, the simplified solution requires 
neither the time of periastron passage to nor the longitude of ascending node, Q (Asada et al. 
2004). If one wishes to know to and Q, however, they can be determined as follows (See also 
Fig. 2). First, we discuss to in this section. 

The projected position of the periastron on the celestial sphere, P^, is determined as 

PA = —{Xe,ye), (14) 

ex 

because the ratio of the semimajor axis to the distance between the center and the focus of 
the ellipse remains unchanged, even after the projection (Asada et al. 2004). The eccentric 
anomaly ua of the periastron in the apparent ellipse is introduced as 

Pa = (acosMyi,6sin'u^), (15) 

where P^ is given also by Eq. (14). Thereby, we can determine (mod 27r). 
By using Eq. (3), we obtain 
S{1,0) _ S{2,1) 

T(1,0) T(2,l)' ^ ^ 

where we can determine 5(1,0) because the eccentric anomaly in the apparent ellipse at to, 
denoted as uo, is nothing but ua, which has been determined by Eqs. (14) and (15). Therefore, 
Eq. (16) is solved for to as 

5(2,0) ^ 5(1,0) ^ 

^° - Wa) ' ~ Wa) " ^ ^ 

where the R.H.S. is obtained from observed quantities. 



4. Longitude of ascending node 



Let us consider the projected periastron at on the apparent eUipse. In the simphfied 
solution, is expressed as Eq. (14). Here we make a translation of {x,y) in such a way that 
the common center of mass can be located at the origin of new coordinates {x',y'). Namely, the 
x' axis is taken to lie along the major axis of the apparent ellipse in the celestial sphere, and 
the y' axis is perpendicular to the x' axis in the celestial sphere (See Fig. 3). In the coordinates 
{x',y'), the position of the projected periastron becomes 

1 — Gi<' 

Pa = -{Xe,ye). (18) 

On the other hand, by projecting the actual ellipse onto the celestial sphere, we obtain 

Pa = (o/<(l — e/^)cosci;,a/<(l — e/^) sintucosz), (19) 

where the coordinates {x,y) are chosen so that the ascending node can be in the x-direction 
(See Fig. 3). 

The longitude of ascending node, which is the angle between the x and x' axes, relates 
the two coordinates of {x',y') and {x,y) by rotation. Therefore, from Eqs. (18) and (19), we 
obtain 

cosil — sinfi \ / axil — gk) cosuj \ 
sinfi cosr2 J \ ax(l — ei^) sincjcosz / 

This relation determines Q (mod 27r). For instance, we obtain explicitly 




(20) 



^ UeCOSU; — XeSmcucosz 
t&nn = — , (21) 

?/e sin cos 2 + Xe cos 

where x^, y^ i and u in the R.H.S. have been determined by Eqs. (5), (6), (8) and (10). 

In conventional methods, determining Q is tightly coupled with u and i. On the other 
hand, it can be done separately from u and i in the simplified solution. 

It should be noted that in practical applications a reference direction chosen by observers 
may be different from the major axis of the apparent ellipse. In such a practical case, Q is the 
angle from the reference direction to the direction of the ascending node. To compute the 
longitude of ascending node, therefore, the angle 6Q from the reference direction to the major 
axis is added into the angle measured from the major axis. In short, the longitude of ascending 
node is generally the sum of Qo and Sfl, where Qo is the angle Q determined by using Eq. (21). 
The expression of SQ is obtained in the straightforward manner, for instance as Eq. (6) in 
Asada et al. (2007), where they denoted Sfl as fl. 

Tables 1 and 2 give an example to show the flow of actual determination of all six 
orbital elements. This would be helpful to the readers who code computing routines in practical 
applications to check their programs. Table 1 shows some given values for all six orbital elements 
and the orbital period. Based on these elements, we flrst prepare a virtual observation data 



5 



set, which is hsted in Table 2. The elements are then reproduced from the data by using the 
proposed method. 

Here, we discuss how to determine in the simplified solution a position of a component 
star at arbitrary time t = t„ (mod T). For ti, t2 and tn, Eq. (3) becomes 
SjnA) ^ 8(2,1) 

T{n,l) T(2,l)- ^ ^ 

This is a transcendental equation for m„ on the celestial sphere. This situation seems similar 
to that the Kepler's equation is transcendental in E on the orbital plane. Here we should note 
that the time of periastron passage is needed in order to treat Kepler's equation, whereas it 
is not for Eq. (22). This is because we employ the time interval T{k,j), while the Kepler's 
equation needs the time itself instead of the interval. Regarding this point, Thiele's method for 
visual binaries is closer to the simplified solution, in the sense that they use the time interval in 
order to delete the time of periastron passage. A crucial difference is that Thiele's method uses 
Kepler's equation on the orbital plane (Thiele 1883), while the simplified one does the constant 
areal velocity in the apparent ellipse on the celestial sphere. In this sense, the simplified solution 
more respects measured quantities on the celestial sphere than conventional ones (See Fig. 2). 
It is verified numerically that the above procedure enables us to determine in the simplified 
solution locations of a star at arbitrary time (See Fig. 4 for an example). 

5. Conclusion 

In this paper, we obtain in the simplified solution both the time of periastron passage and 
the longitude of ascending node in order to complete the solution (See Fig. 2). In conclusion, 
the simplified solution requires neither iterative nor numerical methods when we determine all 
the elements including and Q. It does only when we wish to determine the star's position at 
arbitrary time. 

Before closing this paper, it is worthwhile to mention that Eqs. (17), (21) and (22) can 
be applied to a case of open orbits in the straightforward manner. For open orbits, expressions 
of Xe, ye, ^K, ciK, i and uj have been already derived in the framework of the simplified solution 
(Asada 2007). 

We would like to thank Hiroshi Kinoshita for stimulating conversations. We wish to 
thank Toshio Fukushima for his useful comments on the earlier version of the manuscript. This 
work was supported by a Japanese Grant-in- Aid for Scientific Research from the Ministry of 
Education, No. 19035002. 
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Fig. 1. Schematic figures of tlie actual elliptical orbit (top figure) and the apparent one (bottom one). 
The scmimajor and semiminor axes of the actual elliptical orbit arc denoted by and bk, respectively. 
Those of the apparent one are a and b, respectively. The position of the star at each time is denoted by 
shaded circles (in the top figure) and filled circles (in the bottom figure). The triangles indicate the center 
of mass (COM). In the bottom figure, the projected center of mass is located on the line connecting the 
center of the apparent ellipse and the projected periastron (at to). 
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Actual 

Keplerian Positions on 
Orbit Orbital Plane 



^K,eK,to,T 


— ► 




Projection 


1,0), Q 


/ 








to, Q 

"tk, ^k, Vk 



Measured on 
Celestial Sphere 



Fig. 2. Flow chart of our procedure of orbit determination. The thin arrow denotes purely theoretical 
steps, where we initially assume an actual Keplerian orbit parameterized by {aK-,eK,tQ,T). A star's 
position on the orbital plane at each time tj is projected onto the celestial sphere defined by i, uj and 
n. The thick arrow denotes observational steps, where we start from measuring star's positions on the 
celestial sphere as (tj,Xj). The steps of determining ck, olki T and {i,i-u) have been examined (Asada et 
al. 2004, Asada 2007). The remaining steps of computing to, and (tfc,Xfc) in the simplified solution, 
denoted by the dashed arrow, are discussed in this paper. 
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Periastron 




Actual 

Keplerian 

Orbit 



Fig. 3. Actual Keplerian orbit and apparent ellipse in three-dimensional space. We introduce the incli- 
nation angle i, the argument of periastron lo and the longitude of ascending node 51. These angles relate 
two coordinates (x',y') and {x,y), both of which choose the origin as the common center of mass. Here 
the x' axis is taken to lie along the major axis of the apparent ellipse, while the x-axis is taken to lie along 
the direction of the ascending node. 
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Fig. 4. Example of orbit determination. Here we assume a = 1, e = 0.5 as an apparent ellipse on the 
celestial sphere. Let a star located at m = 60, U2 = 130, U3 = 170, U4 ~ 200 (deg.) at each time ti for 
i = 1,2,3,4. Regarding time, we assume that the time interval T{i + is the same as unity, namely 
ti ^ i — 1 (i.e. tl =0), for simplicity. The simplified solution allows for arbitrary time interval. The 
observed positions of the star are denoted by the filled circle. The quantities determined in the present 
procedure are ck — 0.56, uk — 1-2, i = A2 deg. and uj = 79 deg. We obtain to = 0.04 by Eq. (17). The 
star's position at ^5 — 5 is obtained as W5 = 2.3 x 10^ deg. The location is denoted by the circle. Using 
the determined ck, ax, i and uj, we determine at t^ also by employing the conventional procedure 
(indicated by the thin arrow in Fig. 2). In the latter case, we need to take account of the longitude of 
ascending node, Q.. We obtain = 15 (deg.) by Eq. (21). The results of U5 by both methods agree with 
each other. 
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Table 1. Numerical example of orbital elements. The reference direction is taken along the x-axis, which is different from 
the major axis of the apparent ellipse in this example. Q is the angle measured from the reference direction. 



Gk Bk i UJ to T 

1 0.5 7r/8 7r/9 tt/IO 20 



Table 2. Virtual observational data set based on the orbital elements in Table 1. For simplicity, we assume that observations 
are sampled with the same frequency. 





Xj 




1 


0.372003 


0.838658 


2 


-0.0648698 


0.831542 


3 


-0.404177 


0.696231 


4 


-0.646827 


0.509083 


5 


-0.809209 


0.304280 
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